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A systems analysis has been conducted for a Neptune-Triton Exploration Mission in 
which aerocapture is used to capture a spacecraft at Neptune. Aerocapture uses aerody- 
namic drag instead of propulsion to decelerate from the interplanetary approach trajectory 
to a captured orbit during a single pass through the atmosphere. After capture, propulsion is 
used to move the spacecraft from the initial captured orbit to the desired science orbit. A 
preliminary assessment identified that a spacecraft with a lift to drag ratio of 0.8 was re- 
quired for aerocapture. Performance analyses of the 0.8 L/D vehicle were performed using a 
high fidelity flight simulation within a Monte Carlo executive to determine mission success 
statistics. The simulation was the Program to Optimize Simulated Trajectories (POST) 
modified to include Neptune specific atmospheric and planet models, spacecraft aerody- 
namic characteristics, and interplanetary trajectory models. To these were added autono- 
mous guidance and pseudo flight controller models. The Monte Carlo analyses incorporated 
approach trajectory delivery errors, aerodynamic characteristics uncertainties, and atmos- 
pheric density variations. Monte Carlo analyses were performed for a reference set of uncer- 
tainties and sets of uncertainties modified to produce increased and reduced atmospheric 
variability. For the reference uncertainties, the 0.8 L/D flatbottom ellipsled vehicle achieves 
100% successful capture and has a 99.87 probability of attaining the science orbit with a 360 
m/s AV budget for apoapsis and periapsis adjustment. Monte Carlo analyses were also per- 
formed for a guidance system that modulates both bank angle and angle of attack with the 
reference set of uncertainties. An alpha and bank modulation guidance system reduces the 
99.87 percentile AY 173 m/s (48%) to 187 m/s for the reference set of uncertainties. 
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Nomenclature 


AU 

= 

Astronomical Unit 

C A 

= 

Aerodynamic axial force coefficient 

C n 

= 

Aerodynamic normal force coefficient 

C.G. 

= 

Center of Gravity 

DOF 

= 

degree of freedom 

GRAM 

= 

Global Reference Atmospheric Model 

HYP AS 

= 

Hybrid Predictor-corrector Aerocapture Scheme 

JPL 

= 

Jet Propulsion Laboratory 

L/D 

= 

Lift to drag ratio 

LAURA 

= 

Langley Aerodynamic Upwind Relaxation Algorithm 

POST 

= 

Program to Optimize Simulated Trajectories 

SEP 

= 

Solar electric propulsion 

TPS 

= 

Thermal protection system 

AV 

= 

Velocity addition 

a 

= 

Standard deviation 


I. Background 

A. Neptune-Triton Exploration Reference Mission 

The reference Neptune-Triton exploration mission was designed to provide Cassini and Galileo level explora- 
tion of the Neptune system. 2 The reference mission has a science orbiter to explore the Neptune-Triton system and 
two probes that enter Neptune’s atmosphere 60 degrees apart in latitude. The science orbiter is placed in orbit about 
Neptune using an aerocapture maneuver described below. The science orbit can range between 3896 x 355000 km 
and 3896 x 500000 km. The aerocapture maneuver would be followed by propulsive maneuvers to place the space- 
craft in a phasing orbit such that subsequent maneuvers would establish an orbit that would encounter Triton at regu- 
lar intervals. Triton is then used as a tour engine to vary the orbit’s inclination and line of apsides similar to Cas- 
sini’s use of Titan. The reference orbiter would measure atmospheric, magnetic, and gravity characteristics and per- 
form global imaging of both Triton and Neptune. 

Launch dates studied for the Neptune-Triton Exploration mission range from 2016 to 2019. A February 21, 2017 
launch date was chosen for this study with launch on a Delta IV 4050 Heavy launch vehicle inside a 5m fairing and a 
transit time of 10.25 years 3 . The launch spacecraft configuration consists of the orbiter, two entry probes, and a solar 
electric propulsion module. Five months prior to reaching Neptune, the two probes are released sequentially such 
that both probes’ missions are completed before the orbiter reaches Neptune. 4 Four and a half months prior to reach- 
ing Neptune, a trajectory deflection maneuver is performed to target the entry interface point for aerocapture. Thirty 
minutes prior to entry interface the SEP is jettisoned. After atmospheric entry, the spacecraft executes the aerocap- 
ture maneuver described below to place the orbiter in the exploration orbit and begin a two or more year science 
mission. 

7. Aerocapture Overview 

Aerocapture is a form of aeroassist used to insert a spacecraft into a desired orbit at targets with an atmosphere. 
Aerocapture uses aerodynamic forces to dissipate the hyperbolic approach energy to an energy level needed to reach 
a target apoapsis after making a single pass through the atmosphere. An active guidance system must be used during 
the aeropass to compensate for uncertainties in entry flight path angles, atmospheric density profiles, and aerody- 
namics. After exiting the atmosphere, propulsive maneuvers are required to change the spacecraft’s exit orbital ele- 
ments to that of the desired phasing orbit. These maneuvers include a periapsis raise and any needed adjustments in 
apoapsis, inclination, and longitude of ascending node. The aerocapture maneuver is illustrated in Fig. 1. After the 
proper phasing with Triton is achieved, additional propulsive maneuvers are performed to provide the desired Triton 
encounter strategy. The AV required for the initial periapsis raise and assumed apoapsis for phasing was used as a 
performance metric in this study. 
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II. Simulation of Aerocapture 
Orbit Insertion 


Hyperbolic 
Approach 
/ Trajectory 


A high fidelity 3 DOF simulation 
of the aerocapture maneuver used to 
insert the spacecraft into its phasing 
orbit was developed in the program to 
optimize simulated trajectories, 

POST 9 . The aerocapture trajectory was 
simulated from the navigation delivery 
point, nominally 60 seconds prior to 
atmospheric interface, to atmospheric 
exit. The simulation determined the 
spacecraft’s trajectory through Nep- 
tune’s atmosphere and tracked key de- 
sign parameters such as heat loads, 
deceleration loads, and AV required 
for the periapsis raise and apoapsis ad- 
justments. The simulation was run in a 
Monte Carlo using uncertainties in the delivery point, 
spacecraft aerodynamics, and atmospheric density profiles 
to provide statistical data for the design parameters. 

The simulation incorporated delivered states, aerody- 
namics, guidance, and control models specifically devel- 
oped for the aeroshell and a model of Neptune’s atmos- 
phere as shown in Fig. 2. A brief description of each 
model follows. 

A. Atmosphere 

An engineering type model of Neptune’s atmosphere 
developed at Marshall Space Flight Center provided at- 
mospheric state properties and composition. 6 The model, 
named Neptune-GRAM, is a global reference atmospheric 
model. Neptune-GRAM ’s state and composition properties 
were based on data from the Voyager flyby and stellar oc- 
cultations. 6 Refer to reference 6 for a complete description 
of Neptune-GRAM. 

Neptune’s state property versus altitude relationships 
vary with latitude, season and time of day. These varia- 
tions are represented in Neptune-GRAM by a parameter 
termed Fminmax. Fminmax ranges from -1 to +1 and is 
used to select a state property versus altitude profile for a 
particular latitude, season and time of day. Figure 3 shows 
the maximum range of density versus altitude profiles. 

In the aerocapture simulation, the density versus alti- 
tude relationship was made to vary with latitude by mak- 
ing Fminmax a cosine function of latitude and was made to 
vary with season by adding a constant bias term, Fbias, to 
the latitudinal variation as shown in Eq. 1. 


Atmospheric 
Entry Interface 
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Modulation 

Periapsis 

Apoapsis Adjustment 
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(Science Orbit attained 
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Modulation 


Science Orbit 


Atmosphere Exit 


Figure 1. Illustration of aerocapture maneuvers. 
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Figure 2. Models incorporated into POST simulation 
of Neptune aerocapture. 
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The range of Fbias was set to 0.56 to +0.56 so that Figure 3. Neptune atmosphere model density profiles. 
Fminmax stayed within its -1 to +1 bounds. Figure 4 
shows the variation of Fminmax with latitude. 
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Neptune-GRAM also superimposed high frequency per- 
turbations onto the nominal atmospheric data to represent 
random variations in atmospheric properties. The magnitude 
of the perturbations can be scaled using the Neptune-GRAM 
parameter rpscale. A random perturbation seed value was 
used to generate randomly perturbed density values. Figure 5 
shows randomly perturbed density versus altitude profiles for 
Fminmax = -1,0 and +1. 

B. Aerodynamics 

A high fidelity aerodynamics model of the flat bottom el- 
lipsled aeroshell was developed and incorporated into the 
simulation as an aerodynamic coefficient database. The data- 
base supplied axial and normal force coefficients as a func- 
tion of angle of attack in the hypersonic flight regime. 7 The 
aerodynamic force coefficients were considered constant 
throughout the aeropass. Refer to reference 7 for more infor- 
mation regarding aerodynamics. The aeroshell geometry is 
shown in Fig. 6. 


C. Navigation 

The Neptune-Triton Exploration Mission navigation 
model was provided by JPL. The navigation model deter- 
mined vehicle entry states about a nominal -12.82° entry 
flight path angle and 29.0 km/s entry velocity. The modeled 
navigation system delivered the spacecraft to atmospheric in- 
terface with a 3a dispersion of ±0.51° about the nominal entry 
flight path angle 5 . Refer to reference 5 for a complete descrip- 
tion of the navigation model. Table 1 summarizes the naviga- 
tion data used in the simulation. 

Table 1. Navigation Data 


Entry 

Entry Flight Path 

Entry FPA 

Velocity, km/s 

Angle, deg 

Uncertainty, deg 

29.0 

-12.818 

±0.5108 



Fminmax 

Figure 4. Latitudinal Variation of Fminmax. 
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Figure 5. Sample Neptune-GRAM Monte-Carlo 
density output. 



D. Guidance 

The Hybrid Predictor-corrector Aerocapture Scheme 
(HYP AS) aerocapture guidance algorithm developed at John- 
son Space Center provided autonomous guidance for the 
simulation. 8 The HYP AS algorithm is an analytical control 
algorithm based on drag acceleration and altitude rate error 
for an aeropass through an exponential atmosphere. 

In this study, two attitude control schemes were used by 
HYP AS. In the first, attitude control was limited to bank an- 
gle modulation. Bank angle controlled the rate of de- 
scent/ascent and effected drag through changes in atmos- 
pheric density. In this study, bank angle modulation was de- 
fined as the baseline guidance model. In the second attitude 
control scheme, attitude control included both angle of attack 
and bank angle modulation. Angle of attack modulation was 


(1,404 m) 


5,500 m 


4.096 m 



Figure 6. Aeroshell geometry. 
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used first to modulate the L/D and ballistic coefficient to meet the desired vertical L/D or exit velocity. If changes in 
ballistic coefficient were not sufficient to meet the targets then bank angle modulation was used to further effect 
drag. As the required L/D and drag changed due to high frequency density perturbations, angle of attack modulation 
was used to drive the bank angle to 90°. In this study, the combined bank and angle of attack modulation was de- 
fined as the alpha modulated guidance model. Both the baseline and advanced guidance models used roll reversals 
to maintain the wedge angle between the exit orbit and Triton’s orbit to within 1°. Refer to reference 8 for complete 
descriptions of the bank modulated and alpha modulated guidance models. 

E. Control 

A 3 -DOF Pseudo controller developed at Langley 
Research Center was used to approximate the attitude 
dynamics of a 6 DOF system. The controller analyti- 
cally calculated the time and angular travel required to 
reach the guidance commanded attitude. Once calcu- 
lated, the controller ramped bank angle and/or angle of 
attack to the commanded value at a user specified 
maximum acceleration until the attitude rate reached a 
user specified maximum. The maximum acceleration 
and rates are defined such that the 3 -DOF response is 
a good approximation of the 6-DOF system. This ap- 
proach has provided good agreement with 6-DOF sys- 
tems in previous simulations. Figure 7 shows the bank 
response to a bank command for the 3 -DOF controller. 



Figure 7. Bank Response of 3-DOF Bank Controller 


III. Monte Carlo Analysis 

The vehicle performance was quantified by statistical data from Monte Carlo Analyses. The analyses consisted 
of 2000 individual Neptune aerocapture simulations with random perturbations in arrival states, vehicle aerodynam- 
ics and Neptune’s atmosphere. A Monte Carlo executive script created simulation input files with generated pertur- 
bations and coordinated simultaneous execution of the simulations on multiple processors across multiple com- 
puters. Various post processing scripts were used to determine the statistical parameters for the 2000 simulations in 
each Monte Carlo analysis and to generate plots. 

The Monte Carlo analyses consisted of a reference case and three sensitivity case studies. In the reference case, 
delivery, atmospheric and aerodynamic uncertainties were based on state of the art navigation, current knowledge of 
Neptune atmosphere and computational fluid dynamics analyses respectively. Table 2 lists the uncertainties and dis- 
tribution types used in the Monte Carlo reference case. In the first sensitivity case study, the magnitudes of the high 
frequency random density perturbations were reduced by 50%. In the second sensitivity study, the latitudinal varia- 
tion of Fminmax was removed and the uncertainty in mean density increased. In the third sensitivity study, alpha 
modulation was added to the reference bank modulation approach. The Monte Carlo analyses performed are summa- 
rized in Table 3. 
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Table 2. Monte Carlo Uncertainties 


Category 

Variable 

Nominal 

±3a or min/max 

Distribution 

Delivery State 

X position 

19813.3 km 

From covariance 

Correlated 


Y position 

-16908.2 km 

From covariance 

Correlated 


Z position 

2612.7 km 

From covariance 

Correlated 


X velocity 

-22.953 km/s 

From covariance 

Correlated 


Y velocity 

-13.324 km/s 

From covariance 

Correlated 


Z velocity 

11.316 km/s 

From covariance 

Correlated 

Atmosphere 

Random Pertubation seed 

1 

1 to 9999 

Uniform 


Fbias 

0 

-0.56 to 0.56 

Uniform 

Aerodynamics 

Trim angle of attack 

40.0 

±4.0 

Normal 


C A 

0.349 

±0.048 

Uniform 


C n 

1.771 

±0.120 

Uniform 

Mass Properties 

Axial C.G. (Xcg/L) 

0.51 

±0.50% 

Uniform 


Radial C.G. (Zcg/L) 

-0.0166 

±0.125% 

Uniform 


Table 3. Summary of Monte Carlo Analyses 


Case 

Perturbation Scale 

Fminmax 

Guidance 

Reference 

1.0 

f(latitude) 

Bank angle modulation 

Reduced Density Perturbations 

0.5 

f(latitude) 

Bank angle modulation 

Increased Density Uncertainty 

1.0 

global 

Bank angle modulation 

Alpha Modulation 

1.0 

f(latitude) 

Bank + alpha modulation 


IV. Results 


A. Reference Case 

In the reference case, the reference set of uncertain- 
ties were used in the Monte Carlo analyses. Figure 8 
shows the reference case atmospheric density variation 
in the aerocapture altitudes as a ratio of perturbed den- 
sity to nominal density. The reference case density var- 
ies up to a factor of 2.25. 

The guidance compensates for these variations as 
well as delivery, aerodynamic, and C.G. variations by 
using the spacecraft’s available control authority. Dis- 
persions in apoapsis altitude at atmospheric exit result 
when the spacecraft’s control authority is insufficient to 
compensate for the variations. Figure 9 shows the refer- 
ence case dispersion in apoapsis and periapsis altitude 
and Fig. 10 shows a histogram of exit apoapsis altitude. 

For the reference uncertainties, 100% of the cases 
successfully captured and 66.1% of the cases were 
within the desired science orbit apoapsis bounds. The 
dispersion in apoapsis altitude between the 0.13 percen- 
tile and 99.87 percentile was 461.4E+03 km. The apoap- 
sis percentiles are summarized in Table 4. 



100 1 * * i i 

0 0.5 1.0 1.5 2.0 2.5 3.0 

Density Ratio 


Figure 8. Reference Case Atmospheric Density 
Variation 
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2.0 


Table 4. Apoapsis Percentiles - Reference Case 


Apoapsis Statistics 

Altitude, xlO 3 km 

0.13 percentile 

371.3 

50.00 percentile 

477.9 

99.87 percentile 

832.7 


Impulsive maneuvers were used to adjust the atmos- 
pheric exit orbit to a nominal 3986 x 430000 km orbit. 
The 99.87 percentile AV required to attain the nominal 
orbit and the percentage of cases captured were used as 
performance and robustness metrics. The 99.87 percen- 
tile AV was used rather than a 3a value since the AV dis- 
tribution was skewed. The skewing is a result of AV be- 
ing required for periapsis raise in all cases regardless of 
whether or not the target apoapsis is met. Any error in 
apoapsis only results in increased AV. Table 5 summa- 
rizes these metrics for the reference case. 


Table 5. Performance Metrics - Reference Case 


Robustness Statistics 

% 

Cases Captured 

100 

Cases within Target Bounds 

66.1 

Performance Statistics 

A V, m/s 

0.13 percentile 

88 

50.00 percentile 

141 

99.87 percentile 

360 


The system is sufficiently robust to overcome refer- 
ence delivery, atmospheric and aerodynamic uncertain- 
ties. A 99.87 probability of attaining the science orbit is 
possible with a 360 m/s AV budget for periapsis raise and 
apoapsis adjustment. The reference case AV histogram is 
shown in Fig. 1 1 . 

B. Reduced Atmospheric High Frequency Density 
Perturbations 

In the reduced high frequency density perturbation 
sensitivity study, the magnitude of atmospheric high fre- 
quency perturbations were reduced 50% using the 
GRAM perturbation multiplier rpscale. The reduction 
may be possible with improved knowledge of Neptune’s 
atmosphere. A half scale density perturbation near the 
periapsis of a selected aeropass is shown relative to full 
scale in Fig. 12. Note that a random perturbation is 
equally likely to be below the mean density as above it. 
Figure 13 shows the dispersion in apoapsis and periapsis 
altitude and Fig. 14 shows a histogram of exit apoapsis 
altitude. 
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Figure 9. Apoapsis and Periapsis Altitude Dispersion 
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Figure 12. Comparison of Full and Half Scale High 
Frequency Density Perturbations 


Figure 13. Apoapsis and Periapsis Altitude 
Dispersion - Half Scale High Frequency Density 
Perturbations 


For an atmosphere with half scale high frequency 
density perturbations, 100% of the cases successfully 
captured and 81.8% of the cases were within the apoap- 
sis bounds. In addition, the dispersion in apoapsis alti- 
tude was reduced 221.2E+03 km (47.9%) to 240.2E+03 
km. The apoapsis percentiles are summarized in Table 6. 


Table 6. Apoapsis Percentiles - Half Scale High Fre- 
quency Density Perturbations 


Apoapsis Statistics 

Altitude, xlO 3 km 

0.13 percentile 

412.7 

50.00 percentile 

456.3 

99.87 percentile 

652.9 


The performance statistics are summarized in Table 
7. Reducing high frequency density perturbations 50% 
reduces 99.87 percentile AV to 271 m/s, an 89 m/s 
(24.7%) reduction relative to the reference case. The 
smaller density perturbations can be compensated for 
with less control authority. This allows the guidance to 
improve targeting of the science orbit apoapsis and re- 
duce AV needed for apoapsis adjustment. The AV histo- 
gram is shown in Fig. 15. 


Table 7. Performance Metrics - Half Scale High Fre- 
quency Density Perturbations 


Robustness Statistics 

% 

Cases Captured 

100 

Cases within Target Bounds 

81.8 

Performance Statistics 

A V, m/s 

0.13 percentile 

87 

50.00 percentile 

118 

99.87 percentile 

271 
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C. Increased Mean Density Uncertainty 

In this sensitivity study, the latitudinal variation of 
Fminmax was removed. In addition, the Monte Carlo 
range of Fbias was increased from -0.56 < Fbias <0.56 
to -1 < Fbias <+l. This made Fminmax constant 
throughout a given aeropass and resulted in the largest 
possible uncertainty range of Fminmax, -1 to +1. It also 
produces the maximum uncertainty in mean density from 
one aeropass to another. The maximum range of Fmin- 
max is 33% larger than the Fminmax range of the refer- 
ence case. Figure 16 shows the density ratio in the aero- 
capture altitudes. In the altitude range of 125 km to 300 
km where the majority of hyperbolic approach velocity 
is dissipated, the density variations above the nominal 
are approximately 25% larger than the reference atmos- 
phere while density variations below the nominal are ap- 
proximately 50% larger reductions. 

The larger reduction in atmospheric density impacts 
the guidance’s ability to target the apoapsis. For cases in 
which the density decreases, the spacecraft does not have 
the control authority to overcome the centripedal accel- 
eration and descend into more dense atmospheric re- 
gions. As a result it exits with an apoapsis above the 
bounds. The dispersion in apoapsis and periapsis altitude 
is shown in Fig. 17. Figure 18 shows a histogram of exit 
apoapsis altitude. 

For an atmosphere with increased uncertainties in 
mean density, 100% of the cases successfully captured 
and 76.3% of the cases were within the apoapsis bounds. 
The dispersion in apoapsis altitude increased 498.7E+03 
km (108.1%) to 960.0E+03 km. The apoapsis percentiles 
are summarized in Table 8. 


Table 8. Apoapsis Altitude Percentiles - Increased 
Mean Density Uncertainty 


Apoapsis Statistics 

Altitude, xlO 3 km 

0.13 percentile 

324.9 

50.00 percentile 

448.7 

99.87 percentile 

1284.5 


The system is robust enough to capture 100% of the 
cases with increased density uncertainty. However, the 
99.87 percentile AV increased 96.0 m/s to 456 m/s, a 
26.7% increase relative to the reference case. Table 9 
summarizes the performance statistics. The AV distribu- 
tion is shown in Fig. 19. 
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Figure 16. Atmospheric Density Variation for 
Maximum Range of Fminmax 



Figure 17. Apoapsis and Periapsis Altitude Dispersion 
- Increased Mean Density Uncertainty 
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Figure 18. Apoapsis Altitude Histogram - 
Increased Mean Density Uncertainty 
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Table 9. Performance Metrics - Increased Mean 
Density Uncertainty 


Robustness Statistics 

% 

Cases Captured 

100 

Cases within Target Bounds 

76.3 

Performance Statistics 

A V, m/s 

0.13 percentile 

88 

50.00 percentile 

125 

99.87 percentile 

456 


D. Alpha Modulated Guidance 

In the alpha modulated guidance sensitivity study, a 
guidance algorithm with angle of attack modulation in 
addition to bank modulation was used. The reference 
atmospheric uncertainty assumptions were used with the 
alpha modulated guidance case. The alpha modulated 
guidance improves the spacecraft’ s ability to compensate 
for density perturbations through changes in drag. Angle 
of attack modulation changes the drag more quickly than 
banking thus allowing the spacecraft to respond more 
quickly to density perturbations. The dispersion in 
apoapsis and periapsis altitude is shown in Fig. 20. Fig- 
ure 21 shows a histogram of exit apoapsis altitude. 

For a spacecraft with the alpha modulated guidance 
system, 100% of the cases were successfully captured 
and 98.4% of the cases were placed within the apoapsis 
bounds. The dispersion in apoapsis altitude decreased 
335.5E+03 km (72.7%) to 125.9E+03 km relative to the 
reference case. The apoapsis percentiles are summarized 
in Table 10. 


Table 10. Apoapsis Altitude Percentiles - Alpha 
Modulated Guidance 


Apoapsis Statistics 

Altitude, xlO 3 km 

0.13 percentile 

402.6 

50.00 percentile 

425.1 

99.87 percentile 

528.5 


The alpha modulated guidance performance statistics 
are summarized in Table 11. The alpha modulated guid- 
ance reduced the 99.87 percentile AV by 174.0 m/s 
(48.3%) to 186 m/s relative to the reference case. The 
AV histogram is shown in Fig. 22. 



Figure 19. AV Histogram - Increased Mean Density 
Uncertainty 
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Figure 20. Apoapsis and Periapsis Altitude Dispersion 
- Alpha Modulated Guidance 
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Figure 21. Apoapsis Altitude Histogram - Alpha 
Modulated Guidance 
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Table 11. Performance Metrics - Alpha Modulated 
Guidance 


Robustness Statistics 

% 

Cases Captured 

100.0 

Cases within Target Bounds 

98.4 

Performance Statistics 

A V, m/s 

0.13 percentile 

88 

50.00 percentile 

117 

99.87 percentile 

186 



AV, m/s 

Figure 22. AV Histogram - Alpha Modulated Guidance 


V. Conclusion 

The performance analysis has shown that for the current Neptune atmospheric model, the 0.806 L/D flatbottom 
ellipsled spacecraft is a viable design that captures 100% of the cases and has a 99.87 probability of successfully in- 
serting the orbiter into its science orbit with a 360 m/s AV budget for periapsis raise and apoapsis adjustment. 

The analysis has also shown that the design is viable for atmospheric models with different variability assump- 
tions. For an atmosphere with Fminmax uncertainties 33% larger than that of the reference atmospheric model, the 
spacecraft can attain the science orbit with a 99.87 probability given a 456 m/s AV budget for periapsis and apoapsis 
adjustment. If a better understanding of Neptune’s atmosphere leads to a 50% reduction in high frequency density 
perturbation magnitude, the science obit can be attained with a 99.87 probability given a 271 m/s AV budget for 
periapsis and apoapsis adjustment. The analysis has also shown that with a combined angle of attack and bank angle 
modulation the spacecraft can attain the science orbit with a 99.87 probability given a 186 m/s AV budget for peri- 
apsis and apoapsis adjustment. 
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